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Abstract: We explore the problem of braneworld black holes in the heterotic 
braneworld scenario of Lukas, Ovrut, Stelle and Waldram (LOSW). We show that 
black string solutions are unstable, and demonstrate some unusual asymptotics in 
the linearized metric. We also solve the fully coupled brane and bulk Einstein equa- 
tions, finding an exact, though singular, solution which corresponds to a brane black 
hole in which the branes spike apart at the Schwarzschild radius. 
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1. Introduction 

Large extra dimensions and braneworlds have been an active topic of interest over 
the past decade, with many interesting implications in phenomenology, cosmology, 
and gravity (see PJ for reviews into these various aspects of LED's). While many 
concrete predictions have been based on the explicit models of Arkani-Hamed et. al. 
[@], or Randall Sundrum (RS) [§]; scenarios set in a string theory context such as 
KKLT |Qj or the heterotic braneworld || have also generated new and interesting 
ideas in early universe cosmology. The heterotic braneworld in particular has given 
rise to the ekpyrotic [fj, and cyclic universe picture f7|, which has been the subject 
of some controversy ||. 

The heterotic braneworld is an interesting alternative to type II string theory 
based models, and makes active use of the eleventh dimension to provide the "large" 
extra dimension for our braneworld. It is based on the Horava-Witten M-theory 
compactification ||, in which there is a hierarchical compact ification to 4D with 
the eleventh dimension larger than the remaining 6 spatial dimensions which are 
compactified on a Calabi-Yau (CY) manifold [|l(|. The set-up then mimics the RS 
model, in which the 11th dimension plays the role of the distance normal to the 
brane. The curvature of the CY sources wrapped M5 branes, which in turn warp the 
"fifth" dimension in an analogous fashion to the RS model. 

Despite the apparent similarities between the heterotic model and the RS model, 
the presence of the scalar field in the gravitational sector has a huge impact on the 
strong gravitational properties of the braneworlds. In RS, cosmo logical braneworlds 
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are precisely determinable [11], as the field equations are completely integrable [1^ 



In heterotic M-cosmology however, the presence of the bulk scalar means that explicit 
analytic solutions can only be found by assuming an ansatz for the metric ||13| , and 
there is no "Birkhoff" theorem for the bulk 1 . 

Black holes are the other main test case for strong gravitational solutions, and 
for the RS model, have proved to be rather problematic (see W?\ for a review). 
While the Schwarzschild solution on the brane extends to a black string in the bulk 
this string is unstable |19| , and the exact solution is not known. Furthermore, 



the parallels between the RS model and the gauge/gravity correspondence of string 
theory j20| have led to the idea that a classical bulk solution will correspond to a 



quantum corrected black hole |pl |, although the evidence so far is equivocal |p2]] . For 
the heterotic braneworld though, we do not expect a holographic analogy, indeed, it 
is unclear what sort of black hole solution we can expect. Of course the Schwarzschild 
solution should provide a black string metric - but is this a sensible solution? Here, 
we explore the properties of heterotic brane black holes, determining what features 
the metric should have, and what the constraints on the system are. We show the 
black string is unstable, then calculate the linearized black hole solution. We then 
comment on the full nonperturbative problem, showing how, unlike RS, there is no 
approximate model for a mini black hole metric, and explore a candidate brane plus 
bulk solution with brane spherical symmetry. 



2. Overview of the heterotic braneworld and perturbation 
theory 

In this section, we describe the braneworld setup of Lukas, Ovrut, Stelle and Waldram 
(LOSW) ||, and give the background solution of heterotic M-theory. We then derive 
the linearized Einstein and scalar field equations. 

2.1 Heterotic M-theory 

We use the dimensionally reduced five-dimensional effective action consisting of a 
bulk scalar-tensor gravity, and two boundary branes: 
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where R is the five- dimensional Ricci scalar, g±, is the induced metric on each brane, 



^5 = 87tG 5 is the effective five dimensional Newton's constant, and a is an arbitrary 



1 The "modified Birkhoff theorem" alluded to in jl4| in fact makes a rather restrictive metric 
ansatz, and cannot be taken as a general statement on the bulk spacetime. See |15| ], |16 for general 
analytic analyses of more complex situations. 
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coupling constant, parametrizing the number of units of 4-form flux which thread 
the Calabi-Yau 2 . The boundary branes have equal and opposite tensions and are 
positioned parallel to each other at y = ±y (where y is the transverse direction to 
the brane), and we impose a Z 2 symmetry at the position of each brane. 

The resulting five-dimensional equations of motion following from the action 
(fH|) are 

Gab =2<t>ia 4>,b -^9ab4>, c 4>,c +3g ab a 2 e~ 2(l} 

(2.2) 
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- Qa[5(y + y )-5(y- 2/o)]^W 



□0 = 12a V 2 * - 12a[S(y + y ) - 5{y - y )}^—, (2.3) 
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where Greek indices run over the four braneworld dimensions, and Latin indices run 
over all five spacetime dimensions. 

The heterotic braneworld is given by the solution (writing y = x 11 ) 

= 0(y) = lni7 = ln + (2.4) 

[1 + 6ay \ 

ds 2 = a 2 (y)ri^dx tl dx' / - dy 2 = H^rj^dx^dx" - dy 2 (2.5) 

Note that this solution is in Gaussian Normal (GN) gauge, (g yy = l,g yfl = 0), 
a different gauge from the one originally written down in ||. This is primarily 
for calculational convenience, and a general interacting brane system will need two 
coordinate patches, one for each brane |23]| , in order to correctly encode the physical 



degrees of freedom. We have normalized the warp factor a(y) so that a = 1 on the 
positive tension brane. 

2.2 Perturbations of the braneworld 

Following the usual braneworld procedure (see e.g. |26[), we write the perturbation of 
the metric as g a b —* g a b + h a b, and choose to remain in the GN gauge. In addition, we 
would like to keep the coordinate positions of the branes fixed at ±t/ for calculational 
simplicity, which means that our coordinate system is no longer global, and an explicit 
scalar degree of freedom is introduced into the metric perturbation. The physical 
system consists of the bulk and the two branes, and the physical degrees of freedom 
are therefore bulk fluctuations and (potentially) a degree of freedom corresponding 
to the fluctuation in position of each brane. In the GN gauge, we encode this by 
explicitly performing a gauge transformation 

y-y + ZO^), x^x^+^Lrj^f (2.6) 

4a 



2 Note, for convenience, we are using the conventions of p| for a rather than H. 
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which maintains the GN gauge, but shifts the brane to j/o + /• Since we have two 
branes, we can potentially have two different such gauge transformations, and we 
require two gauge patches, one for each brane [23]. In the overlap, we simply perform 
the relevant shift of the ^-coordinate to match the two patches. Under such a change 
of coordinates, the metric and scalar field change according to their Lie derivatives 
along the coordinate transformation: 

2a 



H = f<t>' 



a 



This gives rise to an explicit scalar component in the perturbation. 

After some algebra, the perturbation equations around the background 
found to be 
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where [D] = [8(y + yo) — S(y — yo)} represents the braneworlds, and have included 
a possible matter perturbation T ^ on each brane. (Note, spacetime indices fi, v... 
are raised and lowered by rj^.) 

The solutions to the homogeneous equations split naturally into a zero mode 
sector, and a massive KK tower of tensor modes. Inspecting ( |2.11|) shows that 
the bulk dependence of the zero mode sector must be either proportional to a 2 , or 



a 4 /2a. The first option of a profile proportional to the warp factor gives 



a transverse-tracefree tensor mode, X/w, as the scalar part of h^ v is a pure 4D gauge 
mode. The other profile is only consistent with a 4D scalar mode, and requires 
the presence of brane fluctuations in order to satisfy the boundary conditions at 
each brane. This gives rise to a situation in which we need two coordinate patches 
p3| , in which we have the GN gauge for each individual brane. The full zero mode 
perturbation is given by 
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which is interpreted as the massless graviton (x), and the radion (/). The ± subscript 
denotes the coordinate patches relevant to each brane. The first coordinate patch 
includes the brane located at y — +y and the second includes the brane located 
at y — —yo- Each coordinate patch is Gaussian Normal with respect to the brane 
it includes (but it does not have to be GN with respect to the other brane). The 
transformation on the overlap is read off from ( |2.13p as: 

/ / 



y -»• y + 



2a 2 



2a 2 



(2.15) 



This represents fluctuations in the interbrane distance, and is a 'breathing mode' for 
the fifth dimension. Note that there is no additional scalar mode coming from the 



i-field, as was originally discussed in |25] 



a = l,yo = 1/12 


18.97 


37.53 


56.16 


a = l,yo = 0.99/6 


10.01 


19 


27.92 


a = l/2,yo = 1/6 


9.58 


18.95 


28.37 


a = l/2,?/o = 0.99/3 


5.02 


9.5 


13.96 


a = 0.1,2/o = 9.9/6 


1.005 


1.9 


2.792 


a = 0.01,2/o = 99/6 


0.1005 


0.19 


0.2792 



Table 1: A table of the first three KK mass eigenvalues for differing values of a and yo- 2/0 
is chosen to be either half, or very nearly its maximal value l/6a. Note that as a — * 0, the 
lowest eigenvalue, tuq oc a, as can be seen from the asymptotic form of the Bessel functions 
in (2.23). By inspection, we see this is well approximated by mg ~ 10a. 



For the massive KK tower, d 2 h iiv = —m 2 h lll/ , and d 2 5(f> = —m 2 5(f), and the scalar 



and tensor modes can be treated separately. Combining ( p79 , 2.10 , 2.12|) in the bulk 



leads to a third order equation for the 'scalar' perturbation: 
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Clearly, 6(f) oc a 6 corresponds to the zero mode already discussed, hence there are 
only two other possible solutions 

(a 6 5(f)) 1 oc J±/${mz) , J-^/^mz) 



where 



(2.17) 
(2.18) 



dy a 5 
a 5a 

is a conformal transverse coordinate. However, ( |2.9| ) and ( 2.12|) imply that at each 
brane 
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While we can balance the coefficients of the fractional Bessel functions to satisfy 
this second equation, the first is not consistent with ( |2.12 ) in a neighborhood of the 
boundary for m 2 ^ 0. 

We are thus left with the massive spin 2 modes hj^ = u m {y)x^w' \ with 
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thus u m oc a 3 J ±4 /5 (mz(y)), and we have 




r 4/5(mz + )J 1/5 (mz) + J_ 4/5 (mz + )J_i /5 (mz)] 
/j 2 4/5 (m2+) + J 4 2 /5 (m2; + ) 



(2.22) 



Note this normalization corresponds to a continuum normalization, and as such is 
only strictly valid in the limit a — > 0, y — > oo, with ay fixed. However, for suf- 
ficiently small a it gives a good working approximation for the computation of the 
propagator, and is far less unwieldy than the exact expression involving hypergeo- 
metric functions. This eigenfunction explicitly satisfies the boundary conditions at 
the '+' brane, however, in order for the boundary conditions at the second ' — ' brane 
to be satisfied, we must have (see also PH): 



(2.23) 



J4/5{mz-)J-4/5{mz + ) = J- 4/5 (mz-)J4/ 5 (mz + ) 



This restricts the values of m allowed in a fashion dependent on a and the background 
interbrane distance 2y$. Some sample mass eigenvalues are shown in Table [l|. 



Figure 1: A plot of the first three KK modes for a = 0.01, we have taken the branes 
(shown in red) to be at ±yo = ±33/2. The masses of the KK modes are given in the final 
line of Table [l|. On the right, for comparison, are the first three eigenfunctions for the RS 
model. 

Like the RS scenario, the behaviour of the massive KK tower is determined by 
Bessel functions, although in this case they are fractional Bessel functions. Figure [I] 
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shows the first three eigenfunctions, and for comparison the corresponding RS profile. 
The heterotic eigenfunctions have a much more regular profile along the y direction, 
which is largely due to the power law, rather than exponential, dependence of the 
warp factor on y. 

Drawing this information together, the Green's function for the spin 2 part of 
the perturbation on the heterotic braneworld is given in the continuum limit by: 

4aa 2 (z)a 2 (z') . , ,,,,,, , N , , N 

3 D (x - x) + dm u m (z)u m (z )D m (x - x ) (2.24) 

< Jo 



G R (x,x') 



This now allows us to compute the effect in the (positive tension) brane of a 
source on the brane 

T+ = S(y-y )T+(xnS^ b (2.25) 

The presence of a source on the brane will in general require the introduction of a 
nonzero scalar perturbation, and once again we need two coordinate patches, one for 
each brane. This is a well known result from RS braneworlds, and is interpreted as 
the brane bending in response to the matter source |26| . 
The general perturbation has the form: 



,6 
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(2.26) 
(2.27) 



where F + is the brane bending term in the coordinate patch of the '+' brane, and 
is given by: 



d 2 F~ 
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(2.281 



Pulling all the information together, we see that the metric in the brane produced 
by matter on the brane is, up to a gauge transformation: 



where 
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and the dilaton is given by: 



t2 ( rn_\ _i_ 72 / m > 



J d A x' D (x - x') 



Thus the brane gravity is a Brans-Dicke theory with u = 1/2. 



(2.29) 



(2.30) 



(2.31) 
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We now explore the possibilities for a heterotic brane black hole. While we 
do not have a complete answer to this problem, there are several approaches using 
both perturbation theory, as well as exact solutions. We start by constructing the 
black string, and exploring its regime of stability. Continuing the theme of linearized 
theory, we compare this with the leading order solution for a point particle on the 
brane. Then we turn to exact approaches, first discussing the possibility of trajecto- 
ries in known bulks before looking at the full axisymmetric problem and presenting 
a possible (singular) solution. 



3. The Black String and Perturbation Theory 



A natural first step in looking for a brane black hole solution is to construct the black 
string: 



7 2 2 

as = a 
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G N M 



dr 2 - r 2 d£l\ 



(3.1) 



However, based on our intuition of cylindrical types of horizon, we expect this black 



string to be unstable | 27|| . The instability of the black string in KK theory occurs 
because there is an unstable massive tensor mode of the 4D Schwarzschild metric, 
therefore, provided a the perturbations of a spacetime allow a separation of the 
perturbation into an effective 4D tensor mode with an orthogonal mass eigenfunction, 
the instability of the string persists to more complicated spacetimes. 

It is not difficult to compute the perturbation equations around the curved 4D 



background, and as in we obtain for the 4D TTF mode: 



(□WV + 2fliV A ')-C + 2^ + g) 2 j 
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(3.2) 



where a(y) takes the appropriate form for the heterotic background ( |2.5p . As this 
is a tensor mode the scalar perturbations are not excited, at least to linear order. 
Thus, we can read off the unstable tensor mode as 



where fajS (t, r] 



is the unstable mode: 



(3.3) 
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and h Q ,hi, h 2 and K are all related via the TTF gauge conditions, and are given in 
The parameter Q depends on the mass m of the longitudinal tensor s-wave, 
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and is found numerically [27], however, for the 4D Schwarzschild metric it is well 
approximated by: 



m 



- m GmM . 



(3.5) 



Clearly there will be an unstable mode if the mass of the black string is (roughly) 
less than 1/2(7^7710, where mo is the minimum eigenvalue permitted for the massive 
tensor tower. This minimum value depends on y and a, but for y ~ l/6a, i.e. if 
the branes are close to their maximal separation, the value is well approximated by 
mo — 10a. Thus the onset of the instability is given by 

1 



G N M < 



20a 



(3.6) 



Once the instability has set in, the 
evolution is similar to the KK string and 
is shown in figure ||. As with the RS 
braneworld, the instability is focussed 
on the brane itself, however, in this case 
the ripples in the event horizon across 
the bulk are more uniform, mirroring 
the behaviour of the transverse eigen- 
functions. This does not give any re- 
liable indication of the nonperturbative Figure 2: The evolution of the black string 
behaviour, however it is reassuring that instability for the nearly marginal case of 
the instability does cluster near the brane, G ^ M = 1 .« = °-°b Va = 99/6. Here only 
rather than having some strong bulk be- 5 eigenfunctions are below the critical mass 

i • t, ,i r i J.1 j value, and the evolution of the event horizon 

naviour. One might therefore expect that 

, , 11111 i , • iiii is exaggerated for clarity. The unperturbed 

the true black hole solution would be lo- 

event horizon is shown as a dotted black line. 

cahzed near the brane. 

One can also use linearized theory 
to obtain a far-field approximation of the black hole metric, by computing the lin- 
earized solution for a point source on the brane 



T, u = MS(v)5(y-y )5y v 

For the scalar, (|2.31| ) immediately gives 

s 2aG 5 M 



(3.7) 
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on the brane. For the tensor, using ( |2. 29| ) , and expanding the Bessel functions in 
u m (0) at small m gives a Newtonian potential of 



v{r) = _]^G ! M( l+ 2'/T[ 



(5ar) 8 / 5 3r[|l 2 



(3.9) 
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This is quite an unusual potential because of the presence of the fractional powers 
of r. Note that it was obtained using the continuum approximation and therefore is 
only really valid for very small values of a, and large values of Uq. 



4. General axisymmetric bulk and the black hole solution 

Turning from perturbation theory, a natural approach towards constructing a brane 
black hole is to use a known bulk and explore possible brane trajectories. Normally, 
one uses a spherically symmetric bulk, taking the brane trajectory through the bulk 
at some nontrivial trajectory 9(r), thereby giving rise to a spherically symmetric 
brane. Typically, although these trajectories exist, they do not correspond to 'empty' 
branes, and energy momentum is required on the brane to source the gravitational 
field. These trajectories are solutions of the brane Tolman-Oppenheimer-Volkoff 
equations (in the context of the RS model, see [28| for work on brane stars and 
TOV, and P| for brane and bulk solutions). 



Unfortunately, with the bulk action given by (|2.1|) , there are no spherically sym- 
metric black hole solutions which are asymptotically flat or (a)dS [[3(J. Indeed, even 



if one asks for only planar symmetric solutions, the system of Einstein equations 



has lost its simplicity, and only limited analytic information can be extracted [|15 
Spherically symmetric geometries exist only for special values of a 2 , and have unusual 



asymptotics pl| . In the case of the heterotic braneworld, the sign of the Liouville 
potential does not permit such a solution, as can be readily seen by attempting to 
solve ( |2.2| ). It may seem strange that there is no black hole solution, since the geom- 
etry has arisen from a compactification from 11 dimensions. However, the wrapped 
5 branes, which give rise to the Liouville potential, mean that the dilaton cannot 
remain fixed in the bulk, and not only remove the possibility of an asymptotically 
flat solution, but imply an anisotropy in any bulk solution. 

It seems therefore that to find a solution, a genuinely axisymmetric bulk metric 
is required. Using coordinate freedom, this metric can be expressed (up to a 2D 



conformal gauge group) as ||T6] 



ds 2 = e Mr > z) dt 2 - e 2x{ - r '^- a ^B{r, z)- 1/2 (dr 2 + dz 2 ) - B{r, z)e- a(r ^ dVi 2 u (4.1) 
This metric has the Einstein equations: 

AB = [2e 2x - Qa 2 Be 2x - 2 ^ a ] B~ 1/2 (4.2) 

A(J + Va ■ ^ = -2a 2 e 2 *- 2 *-°B-V 2 (4.3) 
B 

A X + ^(Va) 2 + i(V0) 2 = - G — 2 2^e 2 ^ 2 ^B^\4A) 

+ l(9 ± a) 2 + \{d^) 2 - 2d ±X ^f = (4.5) 
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and for the dilaton: 

A0 + V0 • ^ = -12a 2 e 2 ^ 2 ^B^ 2 (4.6) 
B 

where A is the 2D Laplacian on (r, z) space, with V as the 2D gradient, and d± = 
d r ± id z . 

This system is similar to the axisymmetric spacetimes explored in [fL6| , where 



the general axisymmetric Einstein equations were derived, then analyzed in detail 
for the case of either spherical symmetry, or a cosmological constant. In each case 
three classes of analytic solution were found. It was noted however, that these were 
specialized solutions, derived assuming some (albeit minimal) metric Ansatz, and 
did not represent the full range of possibilities for the spacetime. 

In this heterotic case, with both spherical symmetry and the bulk scalar field, 
the set of equations is more involved, and like the Einstein axisymmetric problem, 
does not have a general solution generating method. Interestingly however, a simple 
(and commonly used) Ansatz of separation of metric variables gives just one family 
of solutions. Setting 

B = b x (r)b 2 {z) (4.7) 
a = a + <7i(r) + a 2 (z) (4.8) 
X = Xo + Xi(r) +X2(z) (4.9) 
<\> = <l> + Mr) + fa(z) (4-10) 



and inspecting (fO) suggests that e 2x B 1//2 is a function of r, and e 2x 2 ^ a B x l 2 is 



a function of z. Other possibilities are that the roles of r and z are swapped (which 
would result in a rotation of the branes), or both are a function of r (or z). Since it 
is the first option which corresponds to the LOSW vacuum, we will use this in order 
to obtain asymptotically flat braneworld solutions. 

Using the restrictions on \ corresponding to e 2x B~ l l 2 being a function of r, and 
e 2x-20-o-^-i/2 a f unc |;ion of z, the equations of motion give the following expressions: 

B = f(z)/g'(r) (4.11) 

, , In f(z) be „ , 
a = a + ag(r) + — !±± - -( (4.12) 

X = Xo + (2b+ 2 1)a g(r) - i In «/(r) + 1 lnf(z) (4.13) 
<f> = 0o + abg(r) + 2 In /(z) + < (4.14) 

where a, b, and c are arbitrary constants, and 

f dz , , 

W 7 " (4 ' 15) 
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The Einstein equations give a pair of NLDE's for / and g: 
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where a dot denotes d/dz and a prime d/dr. The / equation has the solution 
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^sinh(c/5C)e( fc - 6 )</ 6 
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where f3 2 = |(1 — | — y). The g equation can be integrated by making a change of 
variable: 

J2b+l)ag (4.19) 



which gives 



9{r) 



2E 



:ln 
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p-2E 
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where -E 2 = a 2 (l + b + 55 2 /4), and Vj is of course the standard 4D Schwarzschild 
potential. 

Pulling this information together, we see that the general bulk separable solution 

is: 



ds 2 = fle-f \v s (p)^e- h <dt 2 - V s (p)-^[dp 2 + p(p - 2E)dtf] - f 2 V s (ppdC 2 



e 20 = V s (p)^fe 2 < 



(4.21) 



The LOSW vacuum corresponds to a = c = 0. In this case E = and all the 
nontrivial p-dependence drops out leaving us with 



ds 2 = 2a(r )fll/ dx fl dx u - (2«C) 4 ^C 2 



(4.22) 



setting 6ay = (2a() 3 — 1 recovers the original GN form. 

Taking b = c = 0, but a = £^0 recovers the "uniform black string" solution. 
For b 7^ however, a is no longer equal to E, and the metric and scalar react to the 
"source" , leading to the metric 



ds 2 = a 2 (y) 



2E 
P 



ofl+6) 

2E\ ~ 



dt 2 - ( 1 - —j [dp 2 + p{p - 2E)dtf 



with the scalar given by: 



o 2 



ab 

2E\^ 



I - f) « 12 (y) 



2E\^ , 9 

1 dy 2 

P J 

(4.23) 
(4.24) 
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The y co-ordinate is no longer a GN coordinate because of the variation of g yy in p. 

Turning to a braneworld solution, we introduce branes at ( = (±, and compute 
the extrinsic curvature in order to evaluate the boundary conditions: 

K u = -fl^-^e-* (j± - |) g u (4.25) 
K pp = _/V3 e (6-a K± /3 e -, ^ + ac^ ^ (4 26) 

K ee = -f^e^<^e-t> (k + ^ ^ (4.27) 

Clearly, for a brane solution (energy = tension) we require ac = 0. The Israel 
equations then give 

3a[cosh(c/3C) - P' 1 sinh(c/5C)] = 3a (4.28) 

at either brane. (The sign of the energy term is taken care of by a flip in the sign of 
the extrinsic curvature due to the normal pointing outwards rather than inwards.) 
Obviously if c = this is trivially satisfied, but if instead a = 0, there is only one 
solution to ( |4.28| ), ( = 0, hence it is not possible to have the two brane heterotic 
set-up. 

Therefore we conclude that for the two brane spacetime, we require the bulk 
solution with c = 0, and so the full spacetime is given by (|4.23|) , the scalar field by 
( |4.24| ), and the branes can be set at any fixed ^/-coordinate, which we will once more 
take as y = ±y to compare with the background LOSW vacuum. Restricting to the 
+ brane, the braneworld solution is 

a a(l + b) 

2E\z , 9 / 2E\ 



ds 2 = [ 1 - J) di 2 ~[^- J [ d P 2 + P(P - 2E)dQ 2 ] (4.29) 



ab 



e 2 <t> = [ 1 - ^-j (4.30) 



Note however that the interbrane distance is not a constant: 

"2/0 



D= / dy\g yy \ 1 / 2 = 2y V s (p)^ (4.31) 



yo 

For ab > 0, the interbrane distance decreases as p decreases, eventually closing off 
the extra dimension at p = 2E. For ab < however, the reverse is true, the branes 
move apart until at p = 2E the transverse separation is infinite. 

Although either option yields a legitimate spherically symmetric braneworld so- 
lution, a reasonable approach is to compare this exact solution with the linearized 
result of the previous section: 

= 2 J^l , 4 = (4.32) 
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Expanding ([4.30 ) at large p yields: 



giving: 



ab 
P 

haG^M 



1 - 



2a 



-2/5 



(4.33) 



(4.34) 



Thus, matching to this linearized solution leads to a bulk in which the branes become 
infinitely far apart as the null singularity is approached. 

Thus, allowing for an axisymmetric 
bulk with two branes bounding it, and 
assuming that the metric is separable, 
we have derived the general brane ra- 
dially symmetric solution which asymp- 
totes the LOSW vacuum. Unfortunately 
this solution is singular at p = 2E, how- 
ever, it does look like the Schwarzschild 
solution at large p. The solution very 
much resembles a string solution, how- 
ever, the presence of the Schwarzschild 
potential premultiplying the bulk z co- 
ordinate causes the interbrane distance 
to vary with p, and in fact the "string" 
becomes infinite as it becomes singular. 




5. Discussion 



Figure 3: A plot of the contours of gu for 
the singular black tube solution (a = E = 1). 
Note that far from the 'horizon' the con- 
tours become vertical, i.e. depend on the bulk 
transverse distance, but as p decreases, the 
effect of the source begins to be felt. 



To sum up, we have explored the exis- 
tence of braneworld black holes in the 

heterotic braneworld scenario of Lukas, Ovrut, Stelle and Waldram. We have shown 
how black string solutions are unstable, and that the linearized solution has rather 
unusual asymptotics. Unfortunately it was not possible to construct approximate 
brane stars, as the anisotropic nature of the LOSW vacuum means we have no 
spherically symmetric bulk black hole solutions. However, we were able to construct 
an axisymmetric bulk solution which looks like Schwarzschild at large distances, but 
which is singular as the Schwarzschild radius is approached. 

Interestingly, the lack of a spherically symmetric solution for any value of mass 
means that unlike the RS and ADD models, we are unable to construct even a small 
black hole perturbatively on the brane (such as the solutions considered in ||32|| ) 
which seems somehow paradoxical as one might expect a small black hole to be a 
small perturbation. However, this is really a signal of the different bulk physics. 
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In ADD and RS, the bulk is pure Einstein gravity (with or without a cosmological 
constant) and at smaller scales the brane becomes less and less relevant. In LOSW 
however, even a small perturbation will interact with the scalar field, which is the 
breathing mode of the underlying Calabi-Yau manifold, and thus accesses the higher 
dimensional physics this indicates. 

The existence of these separable axisymmetric solutions is an interesting conse- 
quence of the scalar field in the bulk, for the RS model does not have an equivalent 
solution. The appearance of the Schwarzschild potential to an irrational power is 
reminiscent of the Poincare invariant p-brane solutions in pure gravity P3"] . In that 
case, a string solution in higher dimensions was found, and while the solution was 
dependent on only one variable (the radial distance) the effect of extra dimensions 
was to introduce these irrational powers of the Schwarzschild potential. Here we were 
looking for an (m-symmetric solution, with a warped braneworld interpretation, yet, 
the effect of the extra dimension turns out to be extremely similar. 

The solution which corresponds to the braneworld linearized field at large p has 
the branes diverging as we move in to smaller p. This is an extremely singular config- 
uration with an infinite bulk null singularity. If however, we relax our requirements 
and do not demand agreement with the linearized solution, then we can take b > 0, 
in which case the bulk pinches off at p = 2E, which is perhaps slightly preferable be- 
haviour. Examining the linearized scalar equation, ( |2.31| ), shows that for the branes 
to move together, rather than apart, at the linearized level we require T < 0, in 
other words, M < 3p. This would mean matter with a stiff equation of state, and 
unfortunately does not seem to match the separable solution. 

Obviously the ansatz of separability was a choice, used in order to get an ex- 
act analytic solution, and may be considered to be too restrictive (although it is a 
common ansatz used in finding supergravity solutions). Indeed, looking at the be- 
haviour of the linearized solution across the bulk does not appear to give the same 
dependence as the separable solution, although these are in different gauges. Clearly 
a numerical integration would give a better indication of the true nature of the solu- 
tion. However, in spite of all the unattractive features, it is still interesting that the 
heterotic braneworld does admit an analytic "brane-vacuum" spherically symmetric 
solution, which is the first example of an exact braneworld 'black hole' solution with 
a consistent bulk in 5 dimensions. 



Acknowledgements 

We would like to thank Christos Charmousis for useful discussions. B.M. and A.L. 
acknowledge EPSRC and Durham University fellowships respectively. 



- 15 - 



References 

[1] V. A. Rubakov, Phys. Usp. 44, 871 (2001) [Usp. Fiz. Nauk 171, 913 (2001)] 
[arXiv:hep-ph/0104152]. 

D. Langlois, Prog. Theor. Phys. Suppl. 148, 181 (2003) [arXiv:hep-th/0209261]. 
P. Brax and C. van de Bruck, Class. Quant. Grav. 20, R201 (2003) 
[arXiv:hep-th/0303095]. 

R. Maartens, Living Rev. Rel. 7, 7 (2004) [arXiv:gr-qc/0312059]. 

S. H. Henry Tye, Lect. Notes Phys. 737, 949 (2008) [arXiv:hep-th/0610221]. 

[2] N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998) 
[arXiv:hep-ph/9803315]. 

I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 436, 
257 (1998) [arXiv:hep-ph/9804398]. 

[3] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) 
[arXiv:hep-ph/9905221]. 

L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) 
[arXiv:hep-th/9906064]. 

[4] S. Kachru, R. Kallosh, A. Linde and S. P. Trivedi, Phys. Rev. D 68, 046005 (2003) 
[arXiv:hep-th/0301240]. 

[5] A. Lukas, B. A. Ovrut, K. S. Stelle and D. Waldram, Phys. Rev. D 59, 086001 
(1999) [arXiv:hep-th/9803235]. 

[6] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. D 64, 123522 
(2001) [arXiv:hep-th/0103239]. 

[7] P. J. Steinhardt and N. Turok, Phys. Rev. D 65, 126003 (2002) 
[arXiv:hep-th/0111098]. 

[8] R. Kallosh, L. Kofman and A. D. Linde, Phys. Rev. D 64, 123523 (2001) 
[arXiv:hep-th/0104073]. 

D. H. Lyth, Phys. Lett. B 524, 1 (2002) [arXiv:hep-ph/0106153]. 

R. Brandenberger and F. Finelli, JHEP 0111, 056 (2001) [arXiv:hep-th/0109004]. 

[9] P. Horava and E. Witten, Nucl. Phys. B 460 (1996) 506 [arXiv:hep-th/95 10209]. 
Nucl. Phys. B 475 (1996) 94 [arXiv:hep-th/9603142]. 

[10] E. Witten, Nucl. Phys. B 471, 135 (1996) [arXiv:hep-th/9602070]. 

T. Banks and M. Dine, Nucl. Phys. B 479, 173 (1996) [arXiv:hep-th/9605136]. 
A. Lukas, B. A. Ovrut and D. Waldram, Nucl. Phys. B 532, 43 (1998) 
[arXiv:hep-th/9710208]. 

[11] H. A. Chamblin and H. S. Reall, Nucl. Phys. B 562, 133 (1999) 
[arXiv:hep-th/9903225]. 



- 16 - 



P. Binetruy, C. Deffayet and D. Langlois, Nucl. Phys. B 565, 269 (2000) 
[arXiv:hep-th/9905012]. 

C. Csaki, M. Graesser, C. F. Kolda and J. Terning, Phys. Lett. B 462, 34 (1999) 
[arXiv:hep-ph/9906513]. 

J. M. Cline, C. Grojean and G. Servant, Phys. Rev. Lett. 83, 4245 (1999) 
[arXiv:hep-ph/9906523]. 

P. Kraus, JHEP 9912, Oil (1999) [arXiv:hep-th/9910149]. 

P. Binetruy, C. Deffayet, U. Ellwanger and D. Langlois, Phys. Lett. B 477, 285 

(2000) [arXiv:hep-th/9910219]. 

[12] P. Bowcock, C. Charmousis and R. Gregory, Class. Quant. Grav. 17, 4745 (2000) 
[arXiv:hep-th/0007177]. 

[13] A. Lukas, B. A. Ovrut and D. Waldram, Phys. Rev. D 60, 086001 (1999) 
[arXiv:hep-th /9806022] . 

H. S. Reall, Phys. Rev. D 59, 103506 (1999) [arXiv:hep-th/9809195]. 

A. Lukas, B. A. Ovrut and D. Waldram, Phys. Rev. D 61, 023506 (2000) 

[arXiv:hep-th/9902071]. 

U. Ellwanger, Eur. Phys. J. C 25, 157 (2002) [arXiv:hep-th/0001126]. 
R. L. Arnowitt, J. Dent and B. Dutta, Phys. Rev. D 70, 126001 (2004) 
[arXiv:hep-th/0405050]. 

[14] J. L. Lehners, P. McFadden and N. Turok, Phys. Rev. D 75, 103510 (2007) 
[arXiv:hep-th/0611259]. 
P. McFadden, arXiv:hep-th/0612008. 

J. L. Lehners, P. McFadden and N. Turok, Phys. Rev. D 76, 023501 (2007) 
[arXiv:hep-th/0612026] . 

[15] C. Charmousis, Class. Quant. Grav. 19, 83 (2002) [arXiv:hep-th/0107126]. 

[16] C. Charmousis and R. Gregory, Class. Quant. Grav. 21, 527 (2004) 
[arXiv:gr-qc/0306069]. 

[17] R. Gregory, "Braneworld black holes," arXiv:0804.2595 [hep-th]. 

[18] A. Chamblin, S. W. Hawking and H. S. Reall, Phys. Rev. D 61 (2000) 065007 
[arXiv:hep-th/9909205]. 

[19] R. Gregory, Class. Quant. Grav. 17, L125 (2000) [arXiv:hep-th/0004101]. 

[20] H. L. Verlinde, Nucl. Phys. B 580, 264 (2000) [arXiv:hep-th/9906182]. 
S. S. Gubser, Phys. Rev. D 63, 084017 (2001) [arXiv:hep-th/9912001]. 
E. P. Verlinde and H. L. Verlinde, JHEP 0005, 034 (2000) [arXiv:hep-th/9912018]. 
M. J. Duff and J. T. Liu, Phys. Rev. Lett. 85, 2052 (2000) [Class. Quant. Grav. 18, 
3207 (2001)] [arXiv:hep-th/0003237]. 

[21] T. Tanaka, Prog. Theor. Phys. Suppl. 148, 307 (2003) [arXiv:gr-qc/0203082]. 



-17- 



R. Emparan, A. Fabbri and N. Kaloper, JHEP 0208, 043 (2002) 
[arXiv:hep-th/0206155]. 

[22] A. L. Fitzpatrick, L. Randall, and T. Wiseman, JHEP 0611, 033 (2006) 
hep-th/0608208. 

A. Fabbri and G. P. Procopio, Class. Quant. Grav. 24, 5371 (2007), 0704.3728 
[hep-th]. 

A. Fabbri and G. P. Procopio, "The Holographic Interpretation of Hawking 
Radiation," arXiv:0705.3363 [gr-qc]. 

T. Tanaka, "Implication of Classical Black Hole Evaporation Conjecture to Floating 
Black Holes," arXiv:0709.3674 [gr-qc]. 

L. Grisa and O. Pujolas, JHEP 0806, 059 (2008) [arXiv:0712.2786 [hep-th]]. 
R. Gregory, S. F. Ross and R. Zegers, JHEP 0809, 029 (2008) [arXiv:0802.2037 
[hep-th]]. 

A. Flachi and T. Tanaka, "Vacuum polarization in asymptotically anti-de Sitter 
black hole geometries," arXiv:0803.3125 [hep-th]. 

O. Pujolas, "Strongly Coupled Radiation from Moving Mirrors and Holography in 
the Karch-Randall model," arXiv:0809.0005 [hep-th]. 

H. Yoshino, "On the existence of a static black hole on a brane," arXiv:0812.0465 
[gr-qc]. 

[23] C. Charmousis, R. Gregory and V. A. Rubakov, Phys. Rev. D 62, 067505 (2000) 
[arXiv:hep-th/9912160]. 

[24] O. Seto and H. Kodama, Phys. Rev. D 63, 123506 (2001) [arXiv:hep-th/0012102]. 

[25] P. Brax and A. C. Davis, Phys. Lett. B 497, 289 (2001) [arXiv:hep-th/0011045]. 
G. Cynolter, Mod. Phys. Lett. A 20 (2005) 519 [arXiv:hep-th/0209152]. 
P. Brax, C. van de Bruck, A. C. Davis and C. S. Rhodes, Phys. Lett. B 531 (2002) 
135 [arXiv:hep-th/0201191]. 

[26] J. Garriga and T. Tanaka, Phys. Rev. Lett. 84, 2778 (2000) [arXiv:hep-th/9911055]. 

[27] R. Gregory and R. Laflamme, Phys. Rev. Lett. 70, 2837 (1993) 
[arXiv:hep-th/9301052]. 

R. Gregory and R. Laflamme, Nucl. Phys. B 428, 399 (1994) [arXiv:hep-th/9404071]. 

[28] N. Dadhich, R. Maartens, P. Papadopoulos and V. Rezania, Phys. Lett. B 487, 1 
(2000) [arXiv:hep-th/0003061]. 

C. Germani and R. Maartens, Phys. Rev. D 64, 124010 (2001) 
[arXiv:hep-th/0107011]. 

G. Kofinas, E. Papantonopoulos and I. Pappa, Phys. Rev. D 66, 104014 (2002) 
[arXiv:hep-th/0112019]. 

J. Ovalle, "Non-uniform Braneworld Stars: an Exact Solution," arXiv:0809.3547 
[gr-qc]. 



- 18 - 



[29] S. S. Seahra, Phys. Rev. D 71, 084020 (2005) [arXiv:gr-qc/0501018]. 

C. Galfard, C. Germani and A. Ishibashi, Phys. Rev. D 73, 064014 (2006) 
[arXiv:hep-th/0512001]. 

S. Creek, R. Gregory, P. Kanti and B. Mistry, Class. Quant. Grav. 23, 6633 (2006) 
[arXiv:hep-th/0606006]. 

[30] S. J. Poletti and D. L. Wiltshire, Phys. Rev. D 50, 7260 (1994) [Erratum-ibid. D 52, 
3753 (1995)] [arXiv:gr-qc/9407021]. 

[31] K. C. Chan, J. H. Home and R. B. Mann, Nucl. Phys. B 447, 441 (1995) 
[arXiv:gr-qc/9502042]. 

[32] N. Dadhich, R. Maartens, P. Papadopoulos and V. Rezania, Phys. Lett. B 487, 1 
(2000) [arXiv:hep-th/0003061]. 

R. Casadio and B. Harms, Phys. Lett. B 487, 209 (2000) [arXiv:hep-th/0004004]. 
S. Dimopoulos and G. L. Landsberg, Phys. Rev. Lett. 87, 161602 (2001) 
[arXiv:hep-ph/0106295]. 

P. Kanti, "Reading the number of extra dimensions in the spectrum of Hawking 
radiation," arXiv:hep-ph/0310162. 

[33] R. Gregory, Nucl. Phys. B 467, 159 (1996) [arXiv:hep-th/9510202]. 



- 19 - 



